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Abstract 

We report for the first time a comprehensive study of the fourth rank tensor describing third harmonic 
generation (THG) and electric field induced second harmonic (EFISH) in centrosymmetric material from 
two different viewpoints: Group Theory (GT) and the Simplified Bond Hyperpolarizability Model (SBHM). 
We show that the fourth rank tensor related to THG and direct current (DC) EFISH can be reduced to two 
independent elements whereas SBHM always gives only one, reproducing perfectly well EFISH experi¬ 
mental results in Metal Oxyde Semiconductor (MOS). We argue that it is possible to reduce the fourth rank 
tensor describing EFISH to a third rank tensor and further deliver a classical explanation of EFISH regard¬ 
ing symmetry breaking where the term containing in the potential immediately leads to second harmonic 
generation (SHG). 
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I. INTRODUCTION 


Nonlinear optics has become an intensive research area in physics after the discovery of the 
laser by Townes et al. in 1955 [1] and theoretical study of harmonic generation by several seminal 
papers of Bloembergen et al. [2, 3] in 1962. One of the most intensely studied nonlinear optic 
phenomena is the generation of higher harmonics with frequency 2oj known as second harmonic 
generation (SHG) first investigated by Franken et al. in 1961 [4]. The discovery of third harmonic 
generation soon followed by Maker and Terhune [5] and subsequently by Ward [6]. 

Bloembergen and coworkers reported a study on the phenomena of electric field induced second 
harmonic (EFISH) generation [7] and afterwards several efforts, both theoretically [8] and also 
experimentally [9-11] have been performed to understand the physics of EFISH. Some interesting 
results include the experimental determination of the hyperpolarizabilities for glass [11] and the 
separation of the bulk and surface nonlinear contributions [9]. EFISH is generated by a direct 
current (DC) static field that produces a potential gradient which breaks the initial symmetry of the 
inversion symmetric material. Therefore the nonlinear polarization due to EFISH is characterized 
by a fourth rank susceptibility tensor 

To the best of our knowledge EFISH has never been analyzed from the view point of the Sim¬ 
plified Bond Hyperpolarizability Model (SBHM) and Group Theory (GT) altogether, even though 
the model has been able to successfully predict SHG experimental results of Si surfaces/interfaces 
[12-14]. The SBHM advantage is that the model can describe SHG generated by a Silicon (Si) 
surface with only two fitting parameter [13] which is rather different from the standard tensorial 
approach [15]. In our recent work [16] we have also demonstrated how SBHM and GT are related 
for several crystal orientations such as Si(lll), Si(OOl), and Si(Oll) and here we expand SBHM 
to include THG and EFISH. 

In the merit of this objective, this work is organized as follows: In Section II we discuss group 
theory results concerning the fourth rank tensor that describes third order nonlinear interactions 
particularly THG and EFISH. In Section III we derive the corresponding tensor using the SBHM 
and compare it with GT. We show in Section IV that due to symmetry breaking in Si(OOI) and 
Si(III) facets, the resulting group of symmetry can be explained and interpreted by an effective 
third rank tensor. In Section V, a SBHM simulation is applied to fit the MOS experiment in Ref. 
[10]. Additionally, a classical, ab initio EFISH model is presented in Section VI. Finally, the 
conclusions of this work are briefly stated. 
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II. GROUP THEORY APPROACH 


Contrary to SBHM, which is basically a phenomenological model to calculate the nonlinear 
response of a medium, group theory is a mathematical logical construction that can be applied to 
investigate the symmetry operations allowed by a crystal. Because the physics must be the same 
upon coordinate transformation, the crystal properties, including the (long wavelength) optical 
ones, must remain unchanged under symmetry operations. Due to specific crystal orientation, 
only certain rotations and mirror planes exist. It is well known that rotations and reflections can 
be represented mathematically by matrices. The set of matrices, which contains all the rotations 
and mirror planes allowed for a particular crystal, is called a point group. 

In general, a third-order nonlinear polarization has the form [17]: 


= ^^xifkl{^q,^m,UJn,UJp)Ej{Um)Ek{Un)Ei{Up) 
jkl mnp 


( 1 ) 


where the subindex in the frequencies denote that they can be different. According to group 
theory, the fourth rank tensor xl/ki for bulk silicon in Eq. (1) belongs to a well known point group 
symmetry labelled Oh- The explicit fourth rank tensor for this group of symmetry is similar to the 
tensor given in Refs. [18, 19]: 
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Therefore, a general fourth rank tensor has in total 81 elements. The four rank tensor in Eq. (2), 
consists of a 3 X 3 matrix, where each component of the matrix consist of 3 x 3 matrix elements. 
As will be shown later, not all components are independent [18, 19]. Eor clarity, the notation used 
for representing the fourth rank tensor is explained here. The first index “f” in Xijki corresponds to 
the rows and the second index “j” to the columns in the main matrix (the external one). It follows 
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then that all the elements in the first row and first column of the external 3x3 matrix have snki 
indices, whereas for the second row and third column of the external matrix it will be S 23 ki and so 
on. In the same way the indices “/c” and will correspond to the usual way of labeling a 3 x 3 
matrix, namely the rows and columns in the inner 3x3 matrix respectively. 

As proven in Refs. [18-21] for three possible different frequencies, the fourth rank tensor 
for Oh has at most 4 independent elements but in special physical cases could be reduced to 3 
[18]. However, for a rotating sample, the tensor needs to be correctly rotated using the following 
procedure [18]: 


'^ijkl ^im (0) Fijn (0) ^ko (0) Rip (0) ^mnop 


(3) 


where i?ab(0) is the rotation matrix. One then obtains for a rotation 0 along the 2 ;-axis the following 
tensor in its most general form: 
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0-1 = "S 3223 + "S 3232 + "S 3322 + 3 s 3333 
O2 = -53223 + -53232 + -53322 ~ -53333 

O3 = -53223 + -53232 ~ -53322 ~ -53333 
O4 = -53223 ~ -53232 + -53322 ~ -53333 

O5 = -53223 ~ -53232 — -53322 + -53333 


A. Third Harmonic Generation 

For the case of third harmonic generation (THG) as long as only a monochromatic frequency 
is used the nonlinear polarization in Eq. (1) takes a simpler form: 

P®(3a;) = '^xifki{^^,0J,uj,u)Ej{u)Ek{u)Ei{u) (5) 

jkl 
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where the fields labeled with “j", “/c” and are now undistinguishable. Because the fields are 
the same they cannot be distinguished by the experimenter and the tensor can be "symmetrized". 
Therefore it is allowed to perform intrinsic permutation of the last three indices [17]: 


Xijkl Xijlk Xiklj Xikjl Xilkj Xiljk 


( 6 ) 


When this permutation is performed on Eq. (2) we obtain the following tensor 
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where ci = | (53223 + S3232 + S3322) or S3223 = -S3232 = 53322- Therefore for THG the general 
tensor that previously consists of 4 independent components is now reduced to two indepen¬ 
dent components. 


B. EFISH 

We now consider the case of a static DC field along the 2 ; axis, and assume one monochromatic 
incident field. The nonlinear polarization in Eq. (1) for this particular situation can be stated as: 

Pf\2uj) = ( 8 ) 

jkl 

where in this case has only the 2 : component different from zero. Because of a monochromatic 
incident field the instrinsic permutation can again be applied for the two middle indices “j” and 
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“/c” resulting in the following tensor: 
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where ci = ^ (S 3232 + S 3322 ) Thus the fourth rank tensor describing EFISH consists of 3 
independent elements. One is tempted at this point to apply Kleinman symmetry for the tensor 
in Eq. (9) where for the nonresonant condition the nonlinear susceptibility is independent of the 
frequency permutation and one can freely interchange all the tensor indices [17]. Therefore we 
have S 3232 = <§3223 and we now have only 2 independent components. There are however reports 
that such a treatment for EFISH cannot always be performed [22-24]. 


III. SBHM ANALYSIS OF THG AND EFISH 

SBHM is a classical phenomenological model first constructed by Aspnes and coworkers that 
can be used to calculate the nonlinear response of anharmonic dipoles where, in its most simplified 
version, the electrons are assumed to oscillate only in the direction of the bonds [13]. In the 
case of nonlinear generation inside the bulk, the polarizabilities and hyperpolarizabilities for a 
monochromatic input field have the same value along all bond directions but to surface bonds, a 
different value is assigned due to symmetry breaking in one particular direction. 

However, symmetry breaking in the bulk can also occur if there is a static electric field applied 
in some specific direction. Let us consider the case where a monochromatic electric field is in¬ 
coming on a Si surface that is rotated along the axis. Here, without a DC field, there are no 
SHG contributions from two photon absorption or dipole contribution inside the bulk because the 
Si atomic cell is centrosymmetric (2x4 opposing bonds) and therefore demands to be zero. 
However, there are several other mechanism that can generate second harmonic signals inside the 
bulk [14], one of them is EFISH which breaks this centrosymmetricity and can be evaluated using 
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SBHM. 


In fact, Peng and Aspnes [25] already describe THG using SBHM. They focus their analysis on 
the far field radiation, calculating and discussing in detail the resulting electric field while outlin¬ 
ing the steps followed to get there through SBHM and therefore omitting the information richness 
that can be obtained from the tensorial description of the problem. They also showed that even if a 
transversal contribution from the bonds is considered, this contribution can be expressed as a con¬ 
stant multiplying the electric field plus exactly the same functional form of the contribution along 
the bonds. Thus, for this reason, our analysis can be done only taking account of the contribution 
along the bonds. From the viewpoint of the SBHM the fourth rank tensor takes the form: 

^ uj) ■ h) ® ■ h) ® ( r ^'^ ■ h) ® ( r ^"^ ■ h) ( 10 ) 

^ i=i 

here V is the volume, are the third order hyperpolarizabilities, bj are the unit vectors in the 
direction of the atomic bonds, and is the rotation matrix about the 2 : axis. The summation 
is performed over all eight bonds which are needed to represent the correct response of the con¬ 
ventional cell. If there exist an applied external DC field along a particular direction e.g. in a 
metal oxide semiconductor (MOS) structure [10] then the nonlinear polarization in Eq. (1) can be 
written as: 

P = ■ ■ ■ E(cu) 0 E(cu) 0 E(0) (11) 

where there are two electric fields oscillating at frequency u and one static or DC field. 

Before going further, we are going to show that SBHM always generate a tensor with Kleinman 
symmetry. Indeed it can be shown e.g. using brute force by applying the most general bond vector 
components that the final tensor is always symmetric. This is due to the following fact: the tensor 
is generated by direct product of the same bonds bj 0 bj ® bj ® bj and the hyperpolarizabilities in 
SBHM are just constants. It was already proven by McGilp that SBHM has Kleinman symmetry 
in the particular case of SHG for surfaces [26]. Here, our line of arguing is general and valid for 
any harmonic generation driven with a single driving frequency (SHG, THG, FHG,...) and applies 
for both bulk or surface. Without loss of generality, we are going to take Eq. (10) only with the 
outer product of the vectors, therefore: 

^(3) ^ (g, g, g, ^^2) 

which is a fourth rank tensor and in terms of its components can be expressed as: 

= hqbrbsbt (13) 
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When the indiees g, r, s, and t in the fourth rank tensor in Eq. (13) are permuted, the bond 
eomponents bq,br,bs, and bt are also permuted in the same way but beeause the eomponents are 
just sealars the produet result will always be the same. Therefore any permutation of the subindiees 
generates exaetly the same tensor and this is Kleinman symmetry. In addition, the tensorial 
produet in Eq. (12) was derived from the one dimensional elassieal equation of motion in SBHM 
[13] whieh assumes motion only along the bonds or it was taken for granted that the eleetronie po¬ 
tential is very mueh aligned along the bonds and even if we wish to take aeeount of the transversal 
eontribution, as was mentioned before this is redueed to a eonstant multiplying the eleetrie field 
plus the very same longitudinal eontribution along the bonds [25]. We now proeeed further to 
diseuss the partieular ease of silieon. The diamond bond direetion inside the Si bulk eonsists of 
8 veetors and the fixed system of referenee ean be seen in Eig. L We ehoose for the orientation 
the "standard" eonventional diamond eell beeause the eorresponding material tensor for different 
erystallographie eells in Refs. [18,19] is defined in this way and their representation ehanges when 
the sample is rotated. Therefore we apply the following bond definition [16]: 

“71 2 
-cosf y 
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FIG. 1: Conventional Cell of Silicon 


where 03 is the third order hyperpolarizability for the bulk and the trigonometrie funetions of angle 
were evaluated in order to simplify the expression. Comparison between Eq. (16) and Eq. (4) 
shows similarities in the tensor eomponents. However, even though this proeedure ean in prineiple 
be done for a general angle 0, the result is very eumbersome. Nevertheless, we ean set 0 = 0 and 
eompare it with Eq. (7): 
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It is straightforward to see from Eq. (17) that the SBHM only demands 1 independent 
parameter which is the third order hyperpolarizahility 03. In this way, the undetermined GT 
eonstants from the independent elements ean be determined in terms of physieal values. Beeause 
the tensor in Eq. (17) is general, the bond model thus prediets that THG and EEISH ean also 


be deseribed by only one independent parameter. As we have shown in the previous seetion, 
the group theory fourth rank tensor for THG and EEISH ean be redueed, respeetively, to only 2 
and 3 independent parameters. Therefore it is very possible that the independent tensor elements 
obtained by GT ean be further redueed. It is now up to the experiment or an ab initio theory, to 
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verify if this is generally true. In fact, we will show in Section V, for the case of EFISH with a DC 
field along the 2 : axis how SBHM can fit experimental results using only one hyperpolarizability 
value. 


IV. EFISH REPRESENTED BY A THIRD RANK TENSOR 


In this section, we argue further that SHG due to a DC EFISH in a particular axis can also be 
described using a third rank tensor with an effective susceptibility. The argument is as follows: It 
is known that Oh posseses one of the highest possible point group symmetries in GT. This group 
includes C4, G 3 , C2, Sq and 5*4 axis, as well as, ah and aa mirror planes [18, 21]. Physically, the 
DC field in the direction normal to the [001] plane or along the 2 : axis, breaks the symmetry in 
the conventional cell: the electronic distribution is no longer the same in other directions and all 
the elements of symmetry that transforms in some way the 2 : coordinate are no longer allowed. 
For this reason, only the axis with symmetry group C 4 and two vertical mirror planes remain: the 
resulting point group is € 4 ^. 

This symmetry breaking is very similar to the SHG result obtained for a Si(OOl) surface even- 
though they are not exactly the same. The main difference is that for the Si bulk two neighboring 
tetrahedral elements are required to model the response of the complete conventional cell whereas 
only one tetrahedral is required to describe the surface. Therefore, the point group for this surface 
is C 2 V, as we have discussed in our previous work [16]. Based on this, we can contract the general 
EFISH tensor in Eq. (2) with a unitary vector in the 2 : direction due to the DC field alienation 
along this axis so that we obtain a third rank tensor associated with a Si(OOl) direction: 
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This third rank tensor can be understood as a 9 x 3 matrix divided into three matrices with dimen- 
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sion 3x3. For a general third rank tensor (“i”, “j”, “/c”= 1, 2, 3), the first index “f” corresponds 
to the rows in the main matrix (the external one). Similarly, the indices “j” and “/c” will corre¬ 
spond to the usual way of labeling a 3 x 3 matrix, where the indices ” j” and ” /c” are respectively 
the rows and columns in the inner 3x3 matrix. Moreover, for the third rank tensors and the 
case of SHG, where the two fundamental driving fields are undistinguishable, we can again apply 
intrinsic permutation and therefore S 3322 = S 3232 , which was already sugested by Eq. (9). Fur¬ 
thermore, assuming symmetry in the diagonal of the contracted matrix representation [19] 
we have S 3223 = S 3232 and again both EFISH and SHG third rank tensor now only requires 
two independent parameters. 

We can compare Eq. (18) with the SBHM tensor in Eq. (17) after contracting it with the DC 
field in x direction: 


( 2 ) _ S0‘2eff 

Xijk,SBHM~EFISH ~ 9V 


0 0 1 
0 0 0 
1 0 0 
0 0 0 
0 0 1 
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which gives the same form as Eq. (18). It has to be emphasized that for the case of EEISH the 
third order hyperpolarizability now takes the form Q; 2 e// = a 2 (ix',a;, 0 ) and can generally differ 
from a 2 {uj,u,u). However, if the electron movement can be described by only one resonator, and 
this is very probable below the band gap, the static EEISH polarizibility has to be connected to 
the q; 2 (ci;, lo, 0), just applying an oscillator model. In Section VI we will show how an expression 
for a 2 eff can be obtained from symmetry breaking of the atomic potential. 

Eurthermore, it is also interesting to explore what happens when the DC field is aligned normal 
to the Si(lll) surface because this orientation has been reported experimentally [10] and studied 
theoretically using SBHM [13, 14]. To achieve this, one needs to transform the general tensor to 
another reference frame before doing the contraction. This additional step is necessary because 
the system of reference shown in Eig. 1 is not the same as that used in standard group theory 
textbooks [18, 19], where the x axis corresponds to the higher order of rotation of the point group 
under analysis. To go from the original system of reference, namely in the direction Si(OOl) to the 
direction of Si(l 11), two independent rotations are required. The first is a rotation of vr/d around 
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the z axis and after that a seeond rotation oi (5/2 about the x axis is performed. One then finally 
arrives to the Si(lll) orientation. The transformation is shown in Fig. 2 but for simplieity it is 
shown only aeting on a tetrahedral element; this tetrahedral element ean be seen with the bonds 
loeated at the left eorner of Fig. 1. 



FIG. 2: Pictorial description of a tetrahedral element rotation to go from the Si(OOl) to Si (111) direction, 
a) SBHM configuration, b) after a vr/d rotation around 2 axis and c) a second rotation of /3/2 around x axis 


Matematieally, these transformations ean be applied to the tensor in Eq. (17) using relation 
from Eq. (3). First, a rotation is applied about the z axis and evaluated at 0 = vr/d whieh we label 
(7r/4) and after that a seeond transformation about the x axis for an angle of /?/2 as mentioned 
before labeled R(®)(/3/2). Therefore, after applying these transformations and eontraeted with a 
unitary veetor in the direetion 2 ; (in this ease the e; axis is parallel to the direetion of Si(l 11) as in 
Fig. 2 e), the resulting effeetive third rank tensor for the Si(l 11) orientation is 
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whieh again shows that in SBHM the third rank tensor deseribing EFISH requires only one inde- 
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pendent parameter which is the second order hyperpolarizability. 

Thus in the same way, when the tensor is transformed from the (001) to (111) direction, the 
symmetry elements in Oh are reduced only to the axis parallel with Interestingly from 
the SBHM viewpoint, the symmetry of a Si(lll) orientation is such that even when two tetra¬ 
hedral elements are used to represent the bulk response, the symmetry remains the same as that 
of a Si(lll) surface (represented by only one tetrahedral element) in a similar manner when we 
describe the effective suseptibility in GaAs(lll) in our previous work [16] but this time due to 
EFISH. Therefore, the SBHM model only establishes one additional relation between the inde¬ 
pendent elements in the tensor because the other ones can be derived from the symmetry in the 
crystal. In addition, modeling EFISH only using one independent element in the tensor is not un¬ 
common at all, Kikuchi and Tada using quantum mechanical perturbative calculations [ 8 ] showed 
that under certain conditions it is possible to use only one independent element in the fourth rank 
tensor to describe EFISH. In Section V we show that using both the fourth and third rank tensor 
we can model EFISH perfectly. 


V, SBHM SIMULATION ON MOS EFISH EXPERIMENT 

We now show that SBHM predicts the experimental EFISH result for a pp (first letter polariza¬ 
tion of the in- field, second letter for -out) of Si(l 11) as published by Aktsipetrov and coworkers 
[10] by using only one hyperpolarizability parameter for the fit. This assumption is valid if we 
consider that the SHG is generated mainly by charge separation in the Si-Si 02 interface. Indeed in 
their transmission experiment, the coherence length between fundamental and SH is long or larger 
than the size of the depletion region, therefore in the measured signal phase matching as well as 
the surface effect can be neglected which is in agreement with the arguments of Refs. [18,19]. 
In addition, internal photoemission between such interface has been studied using time dependent 
(TD) SHG as an interface dopant probe [27-30]. 

We apply the following bond orientation for the Si(l 11) bulk [31] which is the same if one starts 
with the diamond unit cell bond orientation in Eq. (12) and perform the necessary transformation 
into a Si(l 11) orientation: 
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where (5 is the angle between each bond and is equal to 2 cos ^ [1/v^] ~ 109.47° for silicon. The 
far field is obtained by [13, 31, 32] 


17i2a; 
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where i is a 3 x 3 unit tensor, is the considered SHG polarization source and k is the unit 
vector of the electric field in the direction of the laboratory observer and has the form 


k 


cos 6 qX + sin OqZ 


(23) 


Following the setup in [10], we set 6^0 = 0 so the electric field is polarized along the x axis 
only. Because the intensity is in arbitrary unit we can assign arbitrary constant for the third order 
hyperpolarizability. The DC bias is varied for every ±0.25V interval around two reference valleys. 
We first simulate the SHG azimuthal feature for a zero external DC using Eq. (17). For a zero bias, 
only the internal field due to the formation of a depletion field contributes to SHG. Without this 
internal field bulk dipole is automatically zero by SBHM simulation because bi = —65, 62 = —be, 
bs = —bj, and 64 = —bg and the radiated fields produced by the bond charges are cancelling out. 

The simulation result for the pp polarization is given in Fig. 3 and shows the expected 6 fold 
pattern at an arbitrary unit SHG intensity of 0.4 (Fig. 3a) where we assume arbitrary constant 
internal DC bias (depletion voltage). The physical explanation according to SBHM is that for a 
normal incidence field the electric field p-polarization is parallel to the x axis thus affecting each 
rotated down-bond 62 , b^, and 64 equally producing a maximal bond radiation when each down 
bond aligns parallel to the x axis (twice when rotated 360°) therefore a sixfold SHG pattern occurs. 
All this occurs naturally using SBHM and one can even look at each single bond contribution 
before summing them up to get a better physical understanding and found that the up bonds bi 
and 65 do not contribute to any SHG intensity for a normal incidence pp field because the bond 
direction is perpendicular to the direction of the electric field. This is in our opinion one of the 
advantages of SBHM because even though limited to the classical view it gives a strong physical 
picture of how SHG is generated by various sources. The result agrees nicely with the experimental 
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FIG. 3: SBHM simulation result for (a) zero external DC bias revealing the sixfold pattern (b) SHG peak 
intensity as a function of the DC bias. The first valley is at 1/ = —0.25V and the second one is at V = 
— 1.5V and is used as the zero volt reference before a forward and backward voltage is applied producing a 
parabolic curve around the two reference voltages. The simulation agrees very well with the experimental 
results in Ref. [10] 

data reported in Fig. 1 of the Ref. [10]. Our simulation further shows that applying an external 
DC bias does not ehange the sixfold pattern but inerease or deerease the SHG intensity. 

We next vary the DC bias for a given referenee voltage. Aeeording to Ref. [10] there are two 
SHG intensity valleys, namely atV = —0.25V with a SHG arbitrary unit intensity of around 0.32 
and V = —1.25V with a SHG intensity of 0.2. Aktsipetrov and eoworkers argue that the absenee 
of a single valley is either due to eleetrie eharge redistribution inside the Si 02 or due to surfaee 
states reeharging . Whatever is the eause, it somehow shifts the zero bias referenee atV = —0.25V 
to C = —1.5V and we adopt this shift in the simulation to produee two symmetrie parabolie fits 
(Fig 3b). Beeause the SHG intensity is obtained by the square of the field two symmetrie parabolie 
eurve around the two valleys are expeeted and is also produeed by the simulation. The simulation 
result eompares well with experiment result in Fig. 2 of Ref. [10]. We will show in Seetion VI, 
how this experiment ean be fully reprodueed using a third rank tensor approaeh. 


VI. CLASSICAL PICTURE OF EFISH IN CENTROSYMMETRIC CRYSTAL 

In this seetion a brief deseription about a elassieal proeess model of EFISH is provided. It is 
well known that the SHG dipole response in bulk Si(lll) without a statie eleetrie field is zero if 
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we neglect the biatomic structure of Si and absorption. Such a system can be modeled as having 
an atomic potential consisting of only even terms [17]. When a static field is directed towards 
the axis it will break the symmetry of the atomic potential. From the view of the hond model 
this static field produces a projection on the three down bonds so that it slightly shifts the 
equilihrium position of the electron away from the former equilihrium position before the 
DC field was applied. Because the lattice position can be assumed to be static in the process this 
change in equilibrium can be seen as breaking the initial symmetry. 


V(r) 



is symmetric relative 
to the lattice 


FIG. 4: Classical picture of EFISH in Si(lll). Left: A static EFISH will shift the equilibrium position 
from the symmetric lattice position and the potential now becomes asymmetric relative to the lattice. Right: 
The projection of the static field on the down bond. 

To the best of our knowledge this picture of the shifting of the potential has never been discussed 
before. The potential of a centrosymmetric material can be expressed as: 

V (r) = muor'^ — mbr^ (24) 

where here we use r as the radial coordinate in SBHM, which is related to ^ as qEdcV sin 7, with 
7 as angle between the upward bond directions onto the z axis. Now a static field along the z axis 
will alter the potential in the form, abbrevating Edc = Edc sin'j 

V (r) = mojQr^ — mbr^ — qE^cr (25) 
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here ujq is the oscillator resonance frequency, q is the electron charge, and b is the parameter that 
characterizes the third harmonic generation (THG) strength. 

This field will produce a shift in the electron potential relative to its initial position, for the 
bond in [111] direction the shift will be larger than for the three backbonds. The force equation is 
then obtained by differentiating the above equation. The new electron equilibrium positions can 
be obtained by assuming br^ « The solution is 


r = Tq 


qEdc 

2 u() 


(26) 


with ro as initial, symmetric position. We now apply a series expansion for the potential V(r) 
around the new electron equilibrium position, r and obtain 


bE},c + 4E: 
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(27) 


(28) 


when viewed in the force equation is proportional to r^, thus: 


'^L//(2^) 


6bEpc 

Uq 


(29) 


It is interesting that based on Eq. (29) the SHG EFISH susceptibility strength depends on the char¬ 
acterizing parameter b of THG and yields immediately zero if Epc = 0 as expected. Increasing 
the DC EFISH voltage will increase the SHG contribution. The linear dependence of Q; 2 e// on the 
DC field is important because it will result in a parabolic feature in the intensity when the EFISH 
field is altered. Furthermore SBHM immediately predicts that under normal incidence only the 
three backbonds contribute to the EFISH signal, despite they experience a rather small shift pro¬ 
portional to sm(19.5), whereas the top bond along the 111 direction does not contribute, because 
it is perpendicular to the polarization of the wave. This argument derived from SBHM predicts that 
for normal incidence on Si(OOl) EFISH should not occur, because the in-plane symmetry V{x, y) 
is not broken, despite the electrons move out along the 2 : direction by the static field. 
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It is thus clear that the expression for the SHG intensity which were applied to plot the exper¬ 
iment using the polarization formula in Seetion V ean also be obtained via the third rank tensor 
in Eq. (20) with an additional rotation freedom about the z axis. Indeed both cases will result in 
a similar intensity formula yielding a 6 fold dependenee for the pp polarization when the applied 
field is nomal to the surfaee: 


Ipp ~ E'jjfja sin^ 30 (30) 

where a ean either be represented by the effeetive first or seeond order hyperpolarizability. There¬ 
fore the experiment in Ref. [10] can be fully refitted to the same aeeuraey using a third rank 
tensor. 

VII. CONCLUSION 

We show that the fourth rank tensor related to THG and DC EFISH nonlinear polarization 
obtained from group theory belongs to the Oh point group and in its most general form eonsists 
of 4 independent parameters which by symmetry arguments can be redueed further to 2 whereas 
SBHM agrees with all the tensor elements different from zero but only requires one independent 
element in the fourth rank tensor which is the third order hyperpolarizability. The fourth rank 
tensor describing DC EFISH in the axis can be further redueed to a third rank tensor, where GT 
shows 2 independent parameters and SBHM requires one. We then demonstrate that the DC MOS 
EFISH experiment in Ref. [10] ean indeed be fitted using only one hyperpolarizability value as 
the fitting parameter. At the end, we presented an important meehanism, a elassieal pieture of the 
symmetry breaking and briefly demonstrate how SHG can arise. 
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